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ABSTRACT 

Context. HD 179949 is an F8V star, orbited by a close-in giant planet with a period of ~ 3 days. Previous studies suggested that the 
planet enhances the magnetic activity of the parent star, producing a chromospheric hot spot which rotates in phase with the planet 
orbit. However, this phenomenon is intermittent since it was observed in several but not all seasons. 

Aims. A long-term monitoring of the magnetic activity of HD 179949 is required to study the amplitude and time scales of star-planet 
interactions. 

Methods. In 2009 we performed a simultaneous optical and X-ray spectroscopic campaign to monitor the magnetic activity of 
HD 179949 during ~5 orbital periods and ~2 stellar rotations. We analyzed the Can H&K lines as a proxy for chromospheric ac- 
tivity, and we studied the X-ray emission in search of flux modulations and to determine basic properties of the coronal plasma. 
Results. A detailed analysis of the flux in the cores of the Can H&K lines and a similar study of the X-ray photometry shows evi- 
dence of source variability, including one flare. The analysis of the the time series of chromospheric data indicates a modulation with 
a ~11 days period, compatible with the stellar rotation period at high latitudes. Instead, the X-ray light curve suggests a signal with a 
period of ~4 days, consistent with the presence of two active regions on opposite hemispheres. 

Conclusions. The observed variability can be explained, most likely, as due to rotational modulation and to intrinsic evolution of 
chromospheric and coronal activity. There is no clear signature related to the orbital motion of the planet, but the possibility that just 
a fraction of the chromospheric and coronal variability is modulated with the orbital period of the planet, or the stellar-planet beat 
period, cannot be excluded. We conclude that any effect due to the presence of the planet is difficult to disentangle. 

Key words. Planet-star interactions. Stars: activity, Stars: magnetic field. 



1. Introduction 

Approximately 20% of the planets detected so far are gi- 
ants (MpSin; > 0.2 My) in tight orbits around their par- 
ent stars (semimajor axis a < 0.1 AU). These planets are 
dubbed Hot Jupiters (HJs). Star-Planet tidal and/or magne- 
tospheric Interaction (SPI) in these systems were claimed 
to be the possible cause of stellar activity enhancements 
dRubenstein & Schaefeill2000t ICuntz et al.ll2000l) . A diagnostic 



of SPI is a signature of stellar activity observed to vary with the 
planet orbital motion, e.g., the Ca II H&K lines, the X-ray emis- 
sion or the optical continuum, which sample the chromospheric, 
coronal and photos pheric levels of the atmosphere, respectively. 
ICuntz et all d2000t) suggested that in the case of tidal interac- 
tions, activity is modulated by half of the orbital period, while in 
the case of magn etospheric interact ions it varies with the orbital 
period. Recently. iFares et aT suggested that the enhance- 
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ment due to magnetospheric SPIs is more likely modulated with 
the beat period of the system, i.e. the synodic period between the 
ste llar rotation and the orbital period. 

IShkolnik et all d2003l I2005L 120081) found the chromospheric 
emission of a few stars hosting HJs to be variable with the planet 
orbital period. Models assuming a magnetic interaction similar 
to that in the Jupiter-Io system fai l to account f or the observed 
chromospheric excess powers (cf. lLanzal 120091 and references 
therein). Other authors suggest that a magnetic field reconfigura- 
tion on a scale comparable with the star-planet separation can be 
triggered by localized reconnection at the boundary betw een the 
stella r corona and the planetary magnetosphere (cf., e.g., lLanzal 
120091 12011b ICohen et al.ll2009l) . Because of the large coronal 
volumes involved in the energy release, these star-planet mag- 
netic interaction (SPMI) models succeed in accounting for the 
observed energy budget. 

Space-based optical photometry has also provided some ev- 
idence of SPI. The variance of CoRoT-2 stellar flux was found 
to be modulated in phase with the plan et orbital period, rea ching 
a maximum near the time of transit dPagano et al.l 12009b . The 
Microvariability and Oscillation of STars (MOST) space tele- 
scope fou nd a spot that persis ted for hundreds of rotation cycles 
on rBoo dWalker et al.H2008l) . while lLanza etail d2009l l20Tlh 
reported hints for SPI in CoRoT-4 and CoRoT-6. 

A statistical survey of the X-ray emission of stars with close- 
in HJs (a < 0.15 AU) suggested that they may be « 2 - 4 times 
more active than stars with more distant pl anets (a > 1.5 AU; 
Kashy ap et al. 2008). On the other hand. [Poppenhaeger et alJ 
d2010l) and iPoppenhaeger & Schmiti d201 11) suggested that the 
coronae of cool stars exhibit a considerable amount of intrinsic 
variability compared to the effects expected from SPI, and that 
the observations of emission enhancements attributed to SPI can 
be explained by selection effect s. Nonetheless, mag netohydro- 
dynamic (MHD) simulations by ICohen et alJ d2009h show that 
SPI may increase the X-ray luminosity. 

Observational studies of HJ hosting stars show that not all 
systems display hints of interaction. Moreover, individual sys- 
tems may show seasons of activity enhancement and seasons 
without such effects. Some systems show a phase lag between 
the su bplanetary longitude and the peak o f the activity enhance- 
ment (IShkolnik et al] 120031 120051 120081) . Different models try 
to account for these phase lags, eithe r considering the c ase 
where the stellar field is a tilted dipole dMcIvor et al. 120061) . or 
adopting a model based on the propagation of Alfven waves 
along the sub-Alfve nic stellar wind flow relative to the planet 
dPreusse et al.l 120061) . or considering a non-pote ntial magneti c 
field co nfiguration for the clos ed coronal field (lLanzal 120081) . 
Finally, ICranmer & Saarl (120071) modeled the Ca 11 H&K light 
curve of a HJ hosting star for different configurations of the stel- 
lar magnetic field based on solar observations, and showed that 
the SPI may be triggered or may not, depending on the specific 
magnetic field geometry. 

HD 179949 (d - 27 pc) is an F8V star, orbited by a giant 
planet with mass Msin; > 0.916 ± 0.076 My, orbital semi- 
major axis a = 0443 ± 0.026 AU, and period of 3.092 days 
dButler et alj|2006l) . The star showed several seasons of activity 
enhancement modulated with the orbital period, with a chromo- 
spheric hot spot leading the planet by ~ 70° - 80°. Seasons with 
no apparent magnetic interaction between th e star and the planet 
dShkolnik et al.fl2008tlGurdemir et 10120121) were also observed. 

We organized a multi-wavelength campaign to observe 
HD 179949 in September and October 2009, and to obtain nearly 
simultaneous optical and X-ray photometry and spectroscopy. 
Our aim was a quantitative characterization of the stellar activ- 



Table 1. Log of ESPaDOnS observations, and S/N per pixel of 
the Ca 11 H&K and the reference Fe 1 /13928.0827 line cores. 



Date at mid exposure 


Orbital 


Exp. time 


S/N H 


S/N K 


S/N FcI 


(MJD-55095) 


phase 


(s) 








4.739 


0.210 


9800 


264 


141 


133 


6.727 


0.854 


9800 


276 


159 


133 


7.732 


0.177 


9800 


242 


141 


1 18 


8.731 


0.500 


9800 


315 


120 


1 10 


9.736 


0.825 


9800 


241 


140 


126 


10.735 


0.148 


9800 


229 


140 


126 


11.750 


0.477 


5040 


108 


59 


54 


12.717 


0.789 


2520 


144 


88 


73 


13.804 


0.141 


5040 


157 


110 


94 


14.731 


0.440 


9800 


174 


100 


90 


16.752 


0.094 


9800 


127 


69 


58 


19.734 


0.059 


9800 


272 


141 


124 



ity to develop models of the star-planet interaction. We remark 
that data series limited to a fraction of the stellar rotation or 
planet orbital period may not provide sufficient information to 
infe r the correct time s cale and origin of the observed phenom- 
ena dMiller et al.ll2012|) : this is not the case of the present study, 
since it is based on optical observations which covered a time 
range of ~ 16 days, i.e. about two stellar rotations and five or- 
bital periods, and X-ray observations spanning more than one 
month. 

Spectropolarimetric data acquire d during our obse rvation 
campaign were already presented bv iFares et al.l (1201 2k In the 
present paper, we focus our attention on the temporal evolution 
of chromospheric and coronal activity, with the analysis of the 
Can H&K lines (AA 3968.487, 3933.673 A) (Sect. 2), and the 
analysis of X-ray data (Sect. [3} obtained in the framework of 
an XMM-Newton Large Program. The main results and conclu- 
sions are reported in Sect. 4. 

2. Optical spectroscopy 

2.1. Observations and data reduction 

2.1.1. ESPaDOnS spectra 

ESPaDOnS is a high-resolution spectropolarimeter installed at 
the 3.6-m Canada-France-Hawaii Telescope (CFHT) on Mauna 
Kea. It provides spectra spanning the whole optical domain 
(3700 to 10000 A) with a resolution of R~65,000 (when used 
in spectropolarimetric mode). We typically collected two spectra 
per night, for a total exposure time of 9800 s (whenever possi- 
ble) on each of 12 nights, between September 24 and October 
9 2009. The d ata were reduced u sing the fully automated tool 
Libre-ESpRIT dDonatiet al.lll997l) . 

In order to derive a common wavelength calibration, we 
choose the highest S/R spectrum as the reference spectrum and 
cross-correlate the other spectra in the series with this reference 
spectrum. When multiple observations per night are available, 
we compute the average spectrum, thus obtaining a master series 
of 12 spectra. The Can H and K spectra taken with ESPaDOnS 
are shown in Fig. [T] and their signal-to-noise ratio is reported in 
the observation log (Table[T]i. 

2.1.2. SARG spectra 

SARG is a cross-dispersed echelle spectrograph installed at the 
Italian 3.6-m Telescopio Nazionale Galileo (TNG) in La Palma, 
Canary Islands. We collected a total of 18 spectra, covering 
the wavelength range 3600-5140 A with a resolution power 
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R=86,000 and an exposure time of 2100 s. Due to weather vari- 
ability, we performed our observations on September 21, from 
October 1 to 3, and from October 6 to 10 2009 (9 nights), col- 
lecting two spectra per night. 

We reduce the spectra using custom IRAF scripts. We per- 
forate the standard steps in extracting the echelle spectra, care- 
fully defining a new set of echelle apertures and new wave- 
length calibrations for each night, thus avoiding the effects of 
instrumental instabilities. We also automatically remove evident 
spikes (due to either cosmic ray hits or hot pixels) in the ex- 
tracted spectra by means of a custom IDL routine. Wavelength 
calibration is performed as discussed in Sect. 12.1.11 

Since the H and K lines cover most of the correspond- 
ing echelle order, no continuum extraction is possible. Thus, 
we normalize these spectra using the continuum extracted from 
ESPaDOnS spectra. 

For each night of observations, we average the two collected 
spectra, obtaining a master series of 9 spectra, unevenly dis- 
tributed between 21st September and 10th October 2009. The 
log of the observations is reported in Table 

2.1.3. du Pont spectra 

We obtained additional high-resolution optical spectra of 
HD 179949 using the Echelle spectrograph mounted on the 2.5- 
m du Pont telescope at Las Campanas Observatory. The instru- 
ment uses a TEK#5 CCD camera with a 2kx2k format and 24 
micron pixels. Simultaneous wavelength coverage extends from 
3700 to 9300 A with a spectral resolution of about 30,000 with 
a l"slit. 



Table 2. Log of SARG observations, and S/N per pixel of the 
Ca ii H&K and of the reference Fe i /13928.0827 line cores. 



Date at mid exposure 
(MJD-55095) 


Orbital 
phase 


Exp. time 

(s) 


S/N H 


S/N K 


S/N Fc , 


11.320 


0.338 


4200 


85 


67 


59 


12.318 


0.660 


4200 


76 


41 


37 


13.318 


0.984 


4200 


42 


44 


39 


16.302 


0.948 


4200 


66 


50 


46 


17.305 


0.273 


4200 


54 


51 


44 


18.332 


0.605 


4200 


55 


43 


39 


19.343 


0.932 


4200 


78 


51 


47 


20.314 


0.246 


4200 


43 


28 


25 



Each stellar exposure is bias-subtracted and flat-fielded for 
pixel-to-pixel sensitivity variations. After optimal extraction, the 
one-dimensional spectra are wavelength calibrated with a Th/Ar 
arc. Finally the spectra are divided by a flat-field response to 
correct for each order's blaze function. 

Our observing run spans 5 consecutive nights from October 
5 to 9 2009, collecting 10 to 50 spectra each night with exposure 
times of ~100 s. We perform wavelength calibration and contin- 
uum normalization as described in Sect. 12.1.11 and Sect. 12.1.21 
respectively. For each night, we compute the average spectrum, 
obtaining a master series of 5 spectra. The log of the observa- 
tions is reported in Table [3] 

2.2. Analysis of the Ca n H&K lines 

In this section, we study the time variability of the chromo- 
spheric activity of HD 179949 through the Can H&K diagnos- 
tics, searching for modulations phased with stellar rotation or the 
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Fig. 2. Deviations from the Can H (left column) and K (right column) average spectrum, for the ESPaDOnS, SARG and du Pont 
series of spectra (from top to bottom, respectively). In each panel, colors code the different nights of observations. The vertical 
dashed lines delimit the core of the lines, as assumed in the present analysis (see the text). 



Table 3. Log of du Pont observations, and S/N per pixel of the 
Can H&K and Fei ,0928.0827 line cores. 



Date 


Orbital 


Exp. time 


S/N H 


S/N K 


S/N M 


(MJD-55095) 


phase 


(s) 








14.506 


0.368 


840 


197 


169 


158 


15.506 


0.691 


5300 


199 


168 


159 


16.506 


0.014 


3800 


134 


98 


86 


17.506 


0.338 


5800 


104 


98 


80 


18.506 


0.661 


3800 


114 


90 


78 



orbital motion of the planet. In the following, we define a proxy 
for the magnetic activity level and we describe the analysis of its 
time modulation. 

2.2.1 . Extraction of the activity proxy 

To look for night-to-night variations in the core profiles due to 
chromospheric activity, first we compute both the H and K aver- 
age spectral profiles for each series of observations (one per in- 
strument), to be taken as representative of the "average star". In 
the computation, we exclude the spectra between MJD=55 111.5 
and 55114 (corresponding to October 6 and 8 2009), because our 



simultaneous X-ray data (see Sect. indicate enhanced flaring 
activity between these dates. 

Then, we compute pixel-by-pixel differences of each nor- 
malized spectrum in the series from the corresponding high- S/N 
average spectrum (Fig. |2}. These differences are generally con- 
sistent with statistical noise, except for the core of the lines, 
where we find larger deviations with respect to the average spec- 
trum. This is a consequence of the fact that the cores of the Ca n 
H&K lines are sensitive to chromospheric activity. Next, we de- 
fine the line core as the spectral range showing the reversal, and 
more precisely as a ~ 1 .5 A window centered on the line center, 
which includes the wavelength range over which we find signif- 
icant deviations from the average spectrum. 

In order to quantify a proxy for stellar activity, we integrate 
the difference spectra over the wavelength ranges defined above. 
The uncertainties on these "Integrated Deviations" (IDs here- 
after) are computed by propagating the noise from the corre- 
sponding spectra. We remark that wider integration domains do 
not provide significantly different integrals, as the residuals out- 
side the cores have null averages; the only effect would be the 
introduction of additional noise in the result. 

We also note that the cross-correlation in wavelength of the 
spectra ensures IDs which are free from Doppler shifts due to the 
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Fig. 3. Correlation between the IDh and IDk indicators. Blue 
diamonds, green triangles and red squares correspond to 
ESPaDOnS, SARG and du Pont measurements, respectively. 
The dashed line is the outlier-resistant best fit of the trend. 
Symbols in gray correspond to spectra collected during three 
days characterized by an enhanced activity level, including a 
prominent X-ray flare (see Sect. |3}, that were not included in 
the correlation analysis. 



orbital motion of the star, and from possible inaccurate wave- 
length solutions due to instrumental misalignment of the spec- 
tra. In particular, if any misalignment among the spectra were 
present, we would find a systematic pattern in the residuals plot- 
ted in Fig.|2]at the positions of the narrow photospheric lines in 
the wings of the H&K lines. Since the residuals away from line 
centers are consistent with randomly distributed noise fluctua- 
tions, this confirms the validity of the cross-correlation described 
m Sect. |2Xj] 

2.2.2. Time series analysis 

The IDs described in the previous section are computed on a nor- 
malized flux scale. This approach allows us to merge the time 
series provided by the three sets of observations into one time 
series per line, consisting of 26 data points. Fig. [3]shows that the 
two proxies, IDh and IDk, are well correlated (Pearson's coeffi- 
cient r > 0.9 and p-values close to 100%). This result allows us 
to sum the two proxies, thus obtaining a merged time series of 
26 data points with a higher S/N, which takes into account the 
contribution of both lines. Six of these data points were obtained 
during a period of enhanced flaring activity, as indicated by the 
X-ray data (Sect. |3J. Since the search for a periodic signal can 
be affected by the presence of isolated flares, we first analyze the 
merged time series by excluding the corresponding data points. 
We also perform the analysis of the full series in Sect. 12.2.31 

As a preliminary step, we compute the Lomb-Scargle peri- 
odogram of the time series (Fig. |4] top panel). We find the high- 
est power at a frequency v\ ~ 0.1 day -1 , corresponding to a 
period of ~ 10 days, and a secondary peak at vi ~ 0.25 day -1 
(~ 4 days). To check whether the secondary frequency V2 is an 
alias of vi, we compute the spectral window of the series (Fig.|4j 
bottom panel). It looks quite noisy with significant sidelobes cor- 
responding to the daily cadence of the observations and its har- 
monics, thus suggesting that aliasing cannot be responsible for 
the secondary peak at 0.25 day -1 . To confirm this prediction, we 
fit our data with a trigonometric polynomial with frequency v\ 
and subtracted it from the series. Then we run the Lomb-Scargle 
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Fig. 4. Top panel - Lomb-Scargle periodogram of the analyzed 
time series. The solid and dashed lines represent the peri- 
odogram for the original and frequency-filtered series, respec- 
tively. Bottom panel - Spectral window of the analyzed time 
series. 

analysis on the residuals, finding that the vi frequency is still 
present after filtering v\ . Thus, V2 is not an alias of v\ . 

To obtain a more robust estimate of the periodicities in the 
time series, we fit the data with a model consisting of a sinu- 
soidal expansion plus a linear trend, the latter representing any 
possible long-term variability: 



ID — q + m t + A„ sin (lo„ t + (/>„) . 



(1) 



where q and m are the parameters of the linear trend, while A,„ 
(jj„ and <p n are the amplitude, angular frequency and initial phase 
of the «-th term of the sinusoidal expansion, respectively. 

We use a non-linea r least-squares estimation method to find 
the best-fit parameters dGallantll 19751) . The advantage of this ap- 
proach over Scargle's analysis is that it takes into account the 
time variation of the data uncertainties, and it provides confi- 
dence intervals for the fitted parameters. Most importantly, it 
performs a fit of different periodic components by adjusting si- 
multaneously their parameters, while in Scargle's approach a 
multiperiodic variability can be detected and characterized only 
iteratively. 

To analyze the time variations of the ID series, we make dif- 
ferent assumptions on the origin of the variability, and we tune 
the model in Eq. ([TJ accordingly. All the assumed models share 
the linear trend and the Fourier expansion of the stellar rotation, 
characterized by the pulsation w rot = 2n/P mt , where P lot is the 
rotation period. They differ in the inclusion of additional sinu- 
soidal terms, chosen to match the various assumptions we wish 
to test. The models considered in our analysis are listed below 
and their properties are summarized in Table 

- Model A: it assumes that the time variations in the ID series 
are due to active regions in the stellar chromosphere, hard- 
rotating with an unknown period P lot ; hence, in Eq. ([JJi the 
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MJD-55095 [day] 

Fig. 5. Best-fits of the data adopting the models listed in Tableland excluding the flaring activity interval. Solid red, green and blue 
lines represent the best fits of the models A, B and C respectively. Symbols and colors of the data points are the same as in Fig. [3] 



(l>„ parameters are the fundamental angular pulsation <jj rot and 
its harmonics; 

- Model B: it assumes that the signal is the sum of two con- 
tributions: the stellar rotation and a possible planet- triggered 
activity enhancement, phased with the orbital motion of the 
planet. In practice, we add to Model A a secondary sinu- 
soid with fixed p eriod P sec = P m \, = 3.092514 days (cf. 
iButler et"alll2006l) : 

- Model C: to investigate the possibility that the SPI signal 
is not exactly synchronized with the orbital motion of the 
planet, we add to Model A a secondary sinusoidal term with 
a free period P sec . 

To evaluate the optimal number of harmonics to be included 
in each model, we adopt the corrected Akaike Information 
Criterior{3 (AICc), which provides a relative measure of the 
goodness of fit of different models when the number of data 
points is comparable with the number of free parameters 
dBurnham & Anderson! 12002b . This criterion not only rewards 
goodness of fit, but it also includes a penalty function that in- 
creases with the number of estimated parameters, thus discour- 
aging overfitting. 

In Table [4] we summarize the output of the AICc: for each 
model we report the optimal number of harmonics «/,, the num- 
ber of fitted parameters «J3> and the best fit periods. The AICc in- 
dicates that Model A has the maximum likelihood, while Model 
B and Model C have confidence levels <1% relative to Model 



1 The AIC is an extension of the canonical Likelihood-Ratio test for 
non-nested models. 

2 For Model A, n p is given by the sum of two parameters from the 
linear trend (q and m), the rotation period P ml , and two parameters for 
each harmonics (amplitude and initial phase), i.e., n p = 2 + 1 + 2 • n; 
for Model B there are two additional parameters: the amplitude and the 
initial phase of the sinusoid with the period fixed at the orbital period of 
the planet; for Model C there are three additional parameters: the am- 
plitude, the period P xc , and the initial phase of the additional sinusoid. 



Table 4. Summary of the best fit models. 



Model" 


n h 


n p 


^rot 


p 

1 sec 


AICc" 




A 


4 


11 


10.8±0.2 






76% 


B 


3 


1 1 


11.2+0.2 




<1% 


43% 


C 


1 


8 


10.9±0.2 


3.6+0.1 


<1% 


13% 


A 


4 


11 


12.0±0.2 




75% 


<1% 


B 


3 


1 1 


12+0.4 




<1% 


<1% 


C 


1 


8 


11±2 


4.2+0.3 




<1% 



<n) The first three lines refer to the censored time series (Sect. \2.2.2l . 

The remaining lines refer to the full time series (Sect. [2723}. 

(6) Relative confidence with respect to the selected model, as given by 

Akaike's criterion. 

(c) Classical^ 2 likelihood. 



A. Moreover, a canonical^ 2 test on Model A gives a confidence 
level of ~75%. The best fits to the time series obtained with the 
different models are plotted in Fig. [5] 

This test significantly rejects models B and C with respect 
to Model A. Specifically, we find that Model B, i.e the model 
with P sec locked at the planetary orbital period, can hardly be 
accepted as a good representation of the data. As a consequence, 
we exclude the scenario where the chromospheric activity is trig- 
gered or enhanced by the interaction with the planet and phased 
with its orbital motion. This is confirmed also by the phase fold- 
ing of our activity proxy with the orbital period of the planet, 
as shown in Fig. [6] we do not find any systematic variation 
in our activity proxy with t he phase of the orbita l motio n of 
the planet, con t rary t o, e.g., IShkolnik etail (|2003[ 120051) and 
iGurdemir et~al1(l2012l) . 

2.2.3. Analysis of the full time series 

In the analysis described so far we have discarded the chromo- 
spheric flux data simultaneous with the X-ray flare that will be 
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Orbital phase 

Fig. 6. ID series vs. the planetary phase. Symbols and colors are 
the same as in Fig. [3] 
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Fig. 7. Best-fits of the the full time series: the solid red, green 
and blue lines represent the best fits of models A, B and C re- 
spectively. Symbols and colors of the data points are the same as 
in Fig. [3] 



discussed in Sect. [3] We now relax this assumption and show 
the results of the analysis performed on the full time series. In 
principle, this is justified by the possibil ity that the fl are be asso- 
ciated with the magnetic SPI (cf., e.g.. lLanzall2012l and refer- 
ences therein). The best-fits with the same models considered in 
Sect. I2.2.21 are shown in Fig. [7] The goodness of the fit is poorer 
than for the censored time series, especially during the flaring 
time interval, and all fits fail to pass the classical x 2 significance 
test. This is probably a consequence of the enhanced short-term 
variability associated with the flare rather than the effect of in- 
creasing the number of data points and corresponding decrease 
of the number of degrees of freedom. 

Nevertheless, if we assume that at least one of the models 
may still provide an acceptable fit to the data, we can apply the 
AICc approach to find the one with the best likelihood. With 
this presumption, the best model turns out to be Model C includ- 
ing the principal and first harmonics of the rotational modula- 
tion. The two periods returned by the fit are P rot =l 1±2 days and 
P sec =4.2+0.3 days. 




15 

55095 (day) 



Fig. 8. ID time series for the Fe i ,13928.0827 line: the derived 
deviations are consistent with the horizontal line at ID=0 and do 
not show any significant time variability. Symbols and colors are 
the same as in Fig. [3] 



The second most likely model is Model A expanded up 
to the second harmonics of P lot , whose best-fit estimate is 
12.0+0.2 days. The relative confidence of Model A with respect 
to Model C is ~75%, i.e. we cannot reject Model A if we accept 
Model C. On the other hand, we can reject Model B, as its confi- 
dence level relative to Model A is close to zero, i.e. we can safely 
exclude the scenario of a SPI phased with the orbital motion of 
the planet. See Table|4]for a summary of the fitted parameters. 

We remark that models A and C provided similar results. As 
a matter of fact, the best-fit Model C turns out with P sec consis- 
tent with the second harmonic of P lot . Moreover, the estimates of 
P rot returned by the two models are consistent within the uncer- 
tainties. In any case, we stress that both models provide poor fits 
to the full time series. This is also indicated by the fact that the 
uncertainties in the fitted parameters are larger than in the case 
of the censored time series considered in Sect. 12.2.21 i.e., with 
the time interval of the X-ray flare removed. 

2.2.4. Significance of the results 

To assess the robustness of our results, first we check whether the 
derived periodicities are due to instrumental effects. As a san- 
ity check, we perform the same analysis on the Fe i /13928.0827 
line, identified in the same echelle order as the Ca n K line. This 
line, being a photospheric feature, is supposed to be weakly af- 
fected by chromospheric activity. Therefore, we do not expect 
any modulation with the rotation of the star or with the orbital 
motion of the planet; if present, any modulation of its IDs is to be 
ascribed to statistical noise or, in the worst case, to instrumental 
or reduction effects. By integrating the difference spectra over 
a 0.3 A width centered on the Fei line, we obtain the merged 
ID series shown in Fig. [8] The data points are consistent with 
a constant line at ID(t)=0 within the uncertainties, i.e. there is 
no evidence of any coherent modulation vs. time. We can thus 
safely exclude that the periodicities found in the Ca n H&K lines 
stem from instrumental effects or flaws in the reduction proce- 
dures. 

To estimate the significance of our results, we also run a sim- 
ple test on our time series: we shuffle the flux and time measure- 
ments to obtain 10,000 fake data samples; then, for each realiza- 
tion, we compute the maximum-likelihood fit for Model A and 
compute the periodicity P rot . This test aims at destroying the true 
periodicities of the series by time-shuffling and, by consequence, 



7 



Scandariato et al.: A coordinated spectroscopic campaign on HD 179949 



Q 
Q. 




Fig. 9. Probability Distribution Function of the period P mt as re- 
covered from 10,000 realizations by time-shuffling, considering 
all the cases (solid black line) or the realizations returning &xt 
lower than 4 times the best-fit x 2 f° r Model A (dot-dashed red 
line; see Sect. 12.2.41 ). The dashed vertical lines bracket the 3<x 
interval around the best-fit period of 10.8 days. 



it returns only frequencies artificially introduced by either ran- 
dom noise or the actual time sampling of the observations. 

In Fig. [9) we plot the derived Probability Distribution 
Function (PDF) of the parameter P lot , as derived by fitting Model 
A and considering two cases: the dot-dashed red PDF has no 
constraint on the final value of the xl of the best fit, thus in- 
cluding also models having a goodness of fit much worse than 
we obtained with Model A; the solid black-lined PDF includes 
only those models that have a xl lower than four times the x 2 
of Model A: this threshold aims at reducing the impact of the 
low-quality best fits on the PDF. 

When we apply the constraint on the x 2 value, the PDF 
shows a significant value for periods longer than ~ 1.0 days, 
with peaks corresponding to an integer or semi-integer number 
of days. The period of 1.0 day corresponds to the first sidelobe 
of the spectral window (cf. Fig. |4] lower panel), while the other 
periodicities correspond to an integer or semi-integer number 
of days, i.e., the subharmonics of the first sidelobe, as a con- 
sequence of the daily cadence of the subsets of observations ob- 
tained from the same observatory. 

To get an estimate of the False Alarm Probability (FAP), we 
integrate the constrained PDF over the 3<x interval centered on 
each best fit value, which corresponds to the p >99% confidence 
interval. The probability that random noise and/or time sampling 
lead to a period consistent within 3cr with the best-fit estimate of 
P I0 t turns out to be FAP~2.2%. Similarly, the probability that a 
periodicity in the range 3.5-4 days randomly arises from the 
noise present in the time series is FAP<3%. 

2.3. Analysis of the Ca n IRT 

The ESPaDOnS spectra also cover the spectral range containing 
the Ca ii Infra-Red Triplet (IRT) located at i8498, 8542, 8662 A. 
One could expect that the low contrast between the variations in 
the line cores due to chromospheric activity and the continuum 
prevents the analysis of any temporal variability. Nonetheless, in 
this spectral region the number of photon counts is remarkably 
higher than the counts in the H&K spectral region (i.e. the blue 
side of the spectrum). Indeed, the S/N in the core of the IRT lines 
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0.000 
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8664 



Fig. 11. Deviations from the average spectrum of the /I8662 A 
line. Colors code the different nights of observations. The verti- 
cal dashed lines delimit the core of the lines, as assumed in the 
present analysis. 



is higher (typically S/N>250) than the S/N of the H&K lines (see 
TablesQ]|2l and©. 

Moreover, variation s in the line profiles o f the IRT have al- 
ready been studied by IShkolnik et al.1 (120081) . finding that the 
IRT-based activity proxies are correlated with the H&K proxies. 
We therefore analyzed the lines in the IRT, and their temporal 
evolution, in the same way as we analyzed the Ca n H&K lines 
inSect.miandl2T2l 

In Fig. [10] we show the time series of the spectra of the 
three components of the triplet (named Tl, T2 and T3, from 
the bluest to the reddest component) as taken by ESPaDOnS, 
while in Fig.[TT]we plot, as an example, the residuals with re- 
spect to the average spectrum for the /18662 A line. We find that 
the residuals integrated over ~0.5 A around the line center are 
not consistent with random noise, and are correlated with the ID 
of the H&K lines defined in Sect. l232l (Fig.[T2l). 

We perform the same analysis as described in Sect. 12.2.21 
We summed the ID series of the three lines, and fit the co-added 
series with Model A (Table [H using the non linear least-squares 
algorithm (Fig. [TBI: the best-fit P rot is fully consistent within the 
error bars with the one obtained in Sect. 12.2.21 

For the sake of comparison, we overplot in Fig. [13]the best 
fit to the H&K data, fixing the frequencies and scaling the am- 
plitudes of the fitted sinusoids so to minimize the sum of the 
squared residuals in the case of the Ca n IRT indicator. The main 
difference between the fitted model and the scaled model is that 
the former more closely follows the data points, while the latter 
gives a smoother fit. This is due to the fact that the fit is based 
on a lower number of data points. As a consequence, the fitted 
model adjusts itself to each single data point rather than to the 
general trend, as it is the case of the scaled model. Nonetheless, 
the fitted model still reproduces the general trend of the scaled 
model, and this is a clear indication that chromospheric variabil- 
ity could be alternatively performed on IRT line measurements 
if provided by ad hoc observational campaigns. 



3. X-ray data 

HD 179949 was observed 15 times with XMM-Newton, be- 
tween Sep 21 and Oct 27, 2009 (see the observation log in 
Table[5]l. The nominal exposure times ranges from about 16 ks to 
68 ks, depending on the actual short-term XMM-Newton sched- 
ule in the A08 observing season. On Sep 21, Oct 3, and Oct 11, 
2009, strictly simultaneous observations with SARG were per- 
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Fig. 10. Spectra of the line cores of the Ca n IRT of HD 179949 taken with ESPaDOns. Different colors code the nights of observa- 
tion. 
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Fig. 12. Correlation between the co-added ID of the H&K lines 
and the IDs of the three lines in the Can IRT (Tl, T2, and T3, 
from bottom to top). In each panel, the dashed line is the outlier- 
resistant best fit of the trend (Pearson's coefficient always >0.9). 
Symbols in gray correspond to the spectra collected during the 
flaring interval indicated by the X-ray observations (Sect. [3), and 
not included in the correlation analysis. 



formed, while in three other occasions (Oct 1, 2, and 6) optical 
spectroscopy with the same instrument was performed within a 
few hours from the X-ray observations. Similarly, simultaneity 
was achieved on Oct 2 and Oct 6 with CFHT, and on Oct 4, 5, 
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Fig. 13. Best-fit (solid line) for the ID proxy of the Can IRT 
with Eq. (Q} Model A, compared with the best-fit of the H&K 
data (dashed line, see text for more details). The gray symbol 
represents the data point not included in the fitting procedure, 
being likely contaminated by the flare activity as discussed in 
Sect. [3] 



and 6 with du Pont. We employed the mos and pn detectors in 
full window mode and with the medium filter in all observations. 

We reprocess all the data with the standard pipeline (SAS 
V9.0), and perform a time screening for high background inter- 
vals (proton flares) that are removed from the following analysis. 
Only four observations are actually affected by such a problem 
(Tabled- We then extract mos and pn source counts from circles 
of 34.5' and 36' radius, respectively, while the local background 
is sampled from a nearby non-overlapping circle of 85' radius 
in all cases. Background-subtracted count rates are reported in 
Table[5]for each observation and instrument. 

We start the X-ray data analysis extracting X-ray light curves 
with a constant bin size of 600 s: in Fig.[l4]we artificially join in 
time all these light curves to give an overall idea of the source 
variability during one month of observing time. The source 
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Table 5. Log of the XMM-Newton observations of HD 179949. 





Start date 


Time" 




MOS 6 




pn 






(2009) 


MJD 


Exp r 


Rate 


Exp r 


Rate 




Obsld 


CUT) 


(d) 


(ks) 


(cnt/ks) 


(ks) 


(cnt/ks) 


Note 


Obsl 


Sep 21 21:08 


55095.99 


18.25 


36.9 ± 1.5 


14.62 


183 ± 4 


Simultaneous TNG obs 


Obs2 


Sep 23 21:17 


55098.07 


30.76 


51.4 ± 1.3 


25.53 


251 ± 3 




Obs3 


Oct 01 14:31 


55105.70 


16.38 


48.7 ± 1.8 


12.97 


228 ± 4 


TNG obs 40 m later 


Obs4 


Oct 02 01:57 


55106.18 


16.37 


54.4 ± 1.9 


7.40 c 


254 ± 6 


simultaneous CFHT obs 


Obs5 


Oct 02 23:46 


55107.07 


13.61 


41.5 ± 1.8 


10.56 


201 ± 4 


TNG obs 4 h earlier; CFHT obs 2 h later 


Obs6 


Oct 03 15:40 


55107.95 


49.85 


43.4 ± 1.0 


42.29 


202 ± 2 


Simultaneous TNG obs; CFHT obs 2 h later 


Obs7 


Oct 04 23:37 


55109.12 


21.85 


c 38.8 ±1.4 


17.74 c 


189 ±3 


Simultaneous du Pont obs; CFHT obs 1 h later 


Obs8 


Oct 05 14:17 


55109.99 


67.10 


38.8 ±0.8 


47.45 c 


189 ±2 


Simultaneous du Pont obs 


Obs9 


Oct 06 23:20 


55111.22 


37.14 


c 52.7 ±1.2 


30.86 c 


255 ±3 


TNG obs 4 h earlier; simultaneous du Pont and CFHT obs 


ObslO 


Oct 11 15:51 


55115.78 


18.94 


47.6 ± 1.6 


15.22 


219±4 


Simultaneous TNG obs, but low-quality data 


Obsll 


Oct 14 22:25 


55119.09 


25.24 


45.8 ± 1.4 


20.72 


214 ±3 




Obsl 2 


Oct 17 07:08 


55121.40 


17.36 


39.5 ± 1.5 


13.83 


198 ±4 




Obsl 3 


Oct 19 09:04 


55123.49 


19.42 


42.6 ± 1.5 


15.64 


212 ±4 




Obsl4 


Oct 22 22:21 


55127.07 


23.26 


48.5 ± 1.5 


19.01 


237 ±4 




Obsl 5 


Oct 27 19:13 


55131.90 


15.39 


47.6 ± 1.8 


12.12 


226 ±4 





Notes. <n) Central time of the observing window. {b) MOS2 data; MOS 1 data are almost identical. w Net exposure time after screening for high 
background levels. 



showed in fact variability by less than a factor 4 overall, on short 
(hours) and long (days-weeks) time scales, including a promi- 
nent flare detected on October 6 2009 (Obs 9). This variabil- 
ity appears to be typical of late-type stars with an intermedi- 
ate coronal activity level, and it is difficult to identify isolated 
flaring events other than in the case of Obs 9. In the same fig- 
ure we marked the times of the strictly simultaneous SARG, 
ESPaDOnS and du Pont observations. 

Then we perform a spectral analysis of each observation, 
in order to evaluate the coronal plasma temperature and the X- 
ray luminosity. The total number of counts in each pn spectrum 
ranges from w 2700 to 7900, and in the range 700-2000 counts 
for each of the mos spectra. With the available photon counting 
statistics, two-component isothermal (2-T) models are adequate 
to describe the source X-ray emission, and to constrain the in- 
terstellar absorption. In practice, we employ the spectra with the 
highest number of total counts (Obs 6, 8, and 9) to check the 
hydrogen column density toward the source, A^h, and the abun- 
dances of some individual elements in the emitting plasma. The 
Nh turns out to be low and compatible with t he value expected 
for ty pical values of the interstellar density dSavage & Mathisl 
119791) and a source at the distance of HD 179949 (d - 27 pc). In 
all the subsequent spectral fitting, we fixed A^h = 10 19 cm -2 . 

Regarding the chemical composition of the coronal plasma, 
we check that non-solar abundances are required in order to ob- 
tain good quality fits of the observed spect ra. In particular, the 
iron abundance is about half solar (using I Anders & Grevessel 
11989 as reference), while t he stellar pho tosphere is metal- 
rich Jvalenti & Fischerl (12005b rFe/Hl=0. 14: fSousa et al.l (l2008h 
Fe/H=[0.21]). The spectra with the best signal-to-noise ratio also 
show that non-solar abundance ratios are required to fit properly 
some prominent line emission complexes: in fact, the best fitting 
models for Obs 6, 8, and 9 suggest O/Fe = 0.3, Ne/Fe = 0.5, 
and Si/Fe = 0.9, while all the other elements have abundances 
essentially equal to the iron one. The iron underabundance by 
a factor of 3 in corona with respect to the photospheric com- 
position is similar t o the well-studied c ase of t Bootis, another 
planet-hosting star (Maggio et al. 2 01 ll) . 



After this preliminary analysis, we decide to fix the abun- 
dance ratios to the value reported above, and to perform 2-T 
spectral fitting with the iron abundance also fixed to 0.5 solafl 
Finally, given the different exposure times of the XMM-Newton 
observations, we split the 4 longest observations (Obs2, Obs6, 
Obs8, and Obs9) in two or three segments so to increase the 
time resolution of the spectral analysis and to ensure a number 
of total counts similar to those in the other observations. This 
choice also ensures similar statistical uncertainties in the best-fit 
parameters. The final results are reported in Table |6] 

The cooler component in the 2-T model has a temperature 
T\ -3-4 MK, while the hotter component has T 2 ~ 5 - 9 MK. 
The ratio of the plasma emission measures of these two compo- 
nents is in the range EM\ /EM2 ~ 1-17. The best-fitting models 
allowed us to compute the broad-band (0.2-2.5 keVfl source X- 
ray luminosity and an average coronal temperature weighted by 
the emission measure. Figure [TBI shows a scatter plot of the X- 
ray luminosity vs. the average coronal temperature. 

3. 1 . Temporal variability of the X-ray luminosity 

We find that the broad-band X-ray luminosity is well correlated 
with the average coronal temperature, a dependence typical of 
solar-type stars, and also considering the same star at different 
observing times. In Fig.[T5lwe show a comparison of HD 179949 
with measurements of coronal temperatu re and luminosity fo r 
the "integrated" Sun at solar maximum (lOrlando et al.l l200ll) . 
and a number of data points for HD 81809, a G-type star moni- 
tored for several years and showing evidence of a solar-like mag- 



3 We verified that the best-fit Fe abundance, if left free to vary, is 
comprised in the narrow range 0.4-0.6 in solar units. 

4 For computing the X-ray luminosity we choose the 0.2-2.5 keV 
band for ease of comparison with similar data taken from the literature. 
In any case, we have verified that an extension of the bandpass to 10 
keV implies a flux increase by less than 1%, for all our best-fitting mod- 
els, because most of the emission is associated with the cooler compo- 
nent, and thermal spectra are characterized by an exponential decrease 
at large photon energies. 
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Fig. 14. X-ray light curves for the 15 XMM-Newton observations (numbered on top), artificially joined in time, with dates in 
September and October 2009 indicated at the bottom. The cumulative exposure time since the start of the first observation is 
reported in abscissa. The upper (black) curve shows the pn count rate, while the lower (red) curve refers to the mosI detector. All 
time bins have the same nominal duration of 600 s. The symbols on top indicate the strictly simultaneous optical spectroscopy 
(triangles for SARG, squares for du Pont, diamonds for ESPaDOnS observations). 



netic cycle both in the c hromospheric (Ca n H&K) and in the 
coronal X-ray emission dFavata et al.l 12008 b. We conclude that 
most of the observed variability of our target is due to chromo- 
spheric and coronal activity, and any effect due to the presence 
of the planet is difficult to disentangle. 

To support this conclusion, in Fig. Q~6] we show the X-ray 
light curve of the data with uncertainties <30%. Following the 
approach introduced in Sect. 12.2.21 and excluding the days with 
flare activity, we find that this time series is fitted with^ , ~2.2 
by considering a linear trend and a sinusoidal modulation with 
P=4. 2+0.1 days. The confidence level for this fit is quite low 
(p -1.6%), while the FAP is ~ 2% (Fig. [16] bottom panel). 

The fitted peri od is close to half the rotation period at inter- 
mediate latitudes dFares et al.ll2012h and can be interpreted with 
the presence of two active regions co-rotating with the star on 
opposite hemispheres. Alternatively, it is consistent with the beat 
period between the stellar equatorial rotation and the planet or- 
bit, but we regard this similarity as a spurious coincidence since 
we do not find any evidence of a SPI signal with the same peri- 
odicity in the chromospheric time series. 

Regarding the flaring event observed on Oct 6 2009, it oc- 
curred within a few hours from simultaneous du Pont and CFHT 



observations which provided optical spectra characterized by en- 
hanced emission in the core of the Can H&K lines. High chro- 
mospheric activity is also indicated by further du Pont and TNG 
observations on Oct 7 and 8 2009. Taken as a whole, these data 
suggest that the X-ray flare occurred during a period of high 
activity on HD 179949, which perdured for ~ 3 days, i.e. one 
full planetary orbit. In the absence of simultaneous X-ray data, 
the enhanced Can signal could have been included blindly in 
the search for a periodic modulation of the chromospheric sig- 
nature, thus introducing significant non-periodic features in the 
analyzed time series. On the other hand, we cannot exclude that 
this event is also related to SPI. 



4. Discussion and conclusions 

We have conducted a nearly simultaneous optical and X-ray ob- 
servational campaign on the planetary host star HD 179949 in 
September-October 2009, aimed at collecting new data for a bet- 
ter understanding of its activity and the possible magnetic star- 
planet interaction. Our observational campaign covers ~5 orbital 
periods and ~2 stellar rotations. The modeling of the residual 
fluxes in the cores of the Ca n H&K lines shows that the most 
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Table 6. X-ray spectral fitting results" . 



Obsld 



Orbital 
phase 



T\ log£M[ 
(MK) (cm" 3 ) 



3 3+o j 

o a+8:t 
i -°-o.i 

2 9 +0 - 3 

3 8^'^ 
3 1 

t 4+8:? 

-0-5 

3.4+"' 
3 4+^ 

-0-5 

3 6 +8 - 3 

t 1+0.3 
—0.3 

^| 

T 7 + 8:3' 

-'■ -2.0 

3 6 +0 - 3 
IT® 

2 9 + " : ' 
3^8 + fJ 

3 3 + <> :2 
3 4 +8i 
9 7+0.5 



?2 

(MK) 



6.2"^ 
8 

6 5+^ 

""'-0.4 

7 9 +T 1 

-13 

6 3 + ° 4 
59^ 
5 



EM, /EM 2 



(10 28 erg/s) 



2 9+rr 

3 0+8: 2 

T 0+0.2 
-Q-2 

3 5 +iii 

—0.2 
^•"-0.2 

3.r'- 9 
3.0 

3 6 m 

3 +0 - 2 
D - -0.2 

3 +0 - 3 
-0.2 

3 +0 - 2 
-'■"-0.2 

o 9+O.2 
J - -0.2 
9 Q+0.5 

^-oa 

o 7+O.3 

o 7+8:3* 
-Q-2 

3 5 +iii 

-Q-2 

3 4 +ii - 2 

llf 

t -3+8:2 

J-J_Q 2 
t 7+0.2 
-'• -0.2 

3 6 +0J 



dof 



Obsl 
Obs2a 
Obs2b 
Obs3 
Obs4 
Obs5 
Obs6a 
Obs6b 
Obs6c 
Obs7 
Obs8a 
Obs8b 
Obs8c 
Obs9a 
Obs9b 
Obs9c 
ObslO 
Obsll 
Obsl 2 
Obsl 3 
Obsl4 
Obsl 5 



0.160 
0.804 
0.862 
0.300 
0.454 
0.743 
0.966 
0.028 
0.088 
0.407 
0.599 
0.674 
0.762 
0.047 
0.093 
0.137 
0.566 
0.627 
0.376 
0.053 
0.209 
0.768 



51.4 
51.6 
51.3 
51.5 
51.4 
51.4 
51.4 
51.4 
51.4 
51.4 
51.2 
51.1 
51.5 
51.5 
51.4 
51.5 
51.4 
51.3 
51.5 
51.4 
51.5 
51.2 



;A:« 



6.3+j? 
6 3 +0 - 6 

"■•'-0.8 

5 5 +0 - 6 
-'-0.3 

5 6 +0 - 6 
J - -0.3 

6 2+ 15 

"■^-4.4 

7.i + ;- 2 
7.0 
7.4: 

5 

6 5 +tl 

—0.8 

6 7 +0 - 9 
' -0.4 
S 4+0.4 



+i!:3 



3.4 

6.7 

1.2 

10.4 

1.5 

3.3 

2.6 

2.4 

3.8 

3.1 

1.2 

0.7 

8.4 

4.1 

1.5 

3.8 

1.5 

1.4 

16.6 

3.2 

3.0 

0.8 



1.4 


55 


1.3 


65 


1.3 


66 


1.6 


57 


1.5 


48 


1.4 


42 


1.1 


58 


0.9 


57 


1.3 


58 


1.2 


69 


1.1 


62 


1.2 


64 


1.6 


61 


1.2 


53 


0.9 


53 


1.5 


53 


1.4 


70 


1.6 


88 


0.9 


57 


1.2 


69 


1.2 


87 


1.8 


56 



Notes. (n) Two-component isothermal model with fixed abundance ratios (O/Fe=0.3, Ne/Fe=0.5, Mg/Fe=0.9, and other element abundances equal 
to Fe = 0.5 Fe Q ). {b) X-ray luminosity in the 0.2-2.5 keV band. 
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Fig. 15. X-ray luminosity vs. coronal temperature for our target, 
HD 179949, and two other comparison stars (see text). Note that 
the temperature was computed as the mean of the two compo- 
nents, weighted by their respective plasma emission measures. 



likely explanation for the observed variability of the chromo- 
spheric flux is the rotational modulation produced by several ac- 
tive regions. The simplest model assumes a modulation with a 
rotation period of ~ 1 1 days and its first and second harmonics. 
Such a period corresponds to the high-latitude rotation period 



of the star as found by iFares et aT I d20ia . who derived its sur- 
face differential rotation from the analysis of Zeeman Doppler 
Imaging maps. The inclination of the stellar rotation axis, as de- 
rived from Doppler Imaging, is « 60°, indicating that one of the 
poles is always in view and suggesting that the amplitude of the 
modulation produced by equatorial active regions is smaller than 
the modulation due to intermediate or high latitude regions, ow- 
ing to the effect of the foreshortening. On the other hand, the 
modulation of the X-ray coronal flux suggests the presence of 
two active regions on opposite hemispheres of the star rotating 
with a period of ~ 9 days, that corresponds to low latitudes. 
These results indicate that the active regions responsible for the 
chromospheric and coronal variabilities are not co-spatial, as of- 
ten observed in the Sun. In HD 179949, the pressure scale height 
of the coronal loops ranges from « 1.5 x 10 10 to « 4.5 x 10 1() cm, 
i.e., it is a significant fraction of the stellar radius (17-52%), im- 
plying that the coronal emission is not necessarily confined to 
levels close to the stellar surface. 

The scenario with a component of the chromospheric and 
coronal variability being phased with the orbital period of the 
planet, or its beat period with respect to the rotation period at 
some stellar latitude, is regarded as unlikely in view of the results 
of our analysis. When we fit the data excluding the time interval 
with enhanced activity (marked by an X-ray flare) and fixing 
the period of the planet, we find a significantly worse fit than 
assuming only rotational modulation (cf. Sect. 12.2.21 . When we 
keep all the data, we find a secondary periodicity at = 4.2 days: 
this value is close to the second harmonic of the stellar rotation 
period P rot ^11 days, but also to the beat period between the 
planet orbital period and the stellar rota tion period, a s expected 
in the case of a magnetic SPI (cf., e.g., Lanza 2012). However, 
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Fig. 16. Top panel - Best fit of the X-ray light curve, using a 
linear term and one sinusoidal modulation. The vertical dotted 
lines bracket the time range likely contaminated by the flare 
activity and excluded in the fitting, as discussed in Sect. 12.2.21 
The fitted period is P =4.2+0.1 days. Bottom panel - Probability 
Distribution Function of the fitted periods recovered after 10,000 
time shufflings. The FAP is ~2%. 



the quality of the fit is poor and it is not possible to discriminate 
between the latter two alternatives. 

The possibility that the intrinsic evolution of the ac- 
tive regions hides the SPI signal cannot be excluded, of 
course. Considering the results in Sect. 12.2.21 the difference 
of ~ 0.5 days between the secondary period of Model C 
CPjec=3.9+0.1 day) and the orbital period of the planet could 
be produced by the intrinsic evolution of the active regions on 
the star or by their shearing owing to the differential rotation. 
Assuming a simple model in which the amplitude of the SPI 
signal is modulated in a sinusoidal fashion with an evolutionary 
period P evo i, we find that the 3.6-day period is obtained when 
P eV oi - 25 days, i.e. slightly longer than the extension of our 
time series and not significantly different from the pole-equator 
lap time of ~ 29 days that com es from the differential rotation 
parameters of lFares et al.l(l2012l) . Unfortunately, it is not possible 
to prove this conjecture because the number of data points is lim- 
ited, so that no significant result can be obtained by splitting the 
time series in two halves and analyzing separately each subseries 



to search for the effects of the intrinsic spot evolution and/or mi- 
gration. We may also note the coincidence between the above 
period differenc e and the d i fferen ce between the polar rotation 
period found bv lFares et al.l (1201 2l) . i.e., 10.3 ± 0.8 days, and the 
rotation period of Model A, i.e., 10.9 + 0.2 days. However, given 
the large uncertainties, these values are compatible with each 
other, and the claim that the active region evolution affects also 
the determination of the rotation period cannot be su pported. 

When comparing our results with those of iFares et all 

one should consider that the present time series is signif- 
icantly longer than the previous one. Based on the smaller data 
set consisting only of the ESPaDOnS spectra. IFares et al" I d2012h 
modeled the chromospheric variability in the Ca n H&K and Ha 
lines by assuming a model identical to our Model C. It provided 
a rotation period Pi = 8.4 days and an additional periodicity 
at Pj = 4.9 days, that was tentatively associated with a beating 
between the orbital motion and the stellar rotation at some inter- 
mediate latitude. However, the present result, based on a longer 
time series and the removal of the data points affected by flaring, 
points toward a modulation produced by stellar rotation without 
any evidence of a planet-induced component. Since the chro- 
mospheric signatures of the star-pla net interaction are pr esent 
in approximately 50% of the seasons (IShkolnik et al.ll2008l) . our 
apparently negative result is in line with the expected frequency. 
Further observations would be valuable to confirm the presence 
of a signature moving in phase with the orbital motion of the 
planet as f ound in 2001-2002 and 2005-2006 bv lShkolnik etall 
d2008l) and lGurdemir et afl (|2012|) . The transition from an on to 
an off state of chromospheric SPI could be associated with the 
variation of the stellar field configuration along an activity cy- 
cle. A tentative period of a 4-5 years could be derived from 
the above results. Howev er, the lack of an observed signature 
in 2009, also confirmed bv lPoppen haeger et ail (1201 lb . indicates 
that the on-off cycles may not be periodic. As a matter of fact, the 
numerical simulations by ICohen et al.l (1201 Ol 1201 ll) show that 
the appearance of a signature in the stellar atmosphere is asso- 
ciated with the planet's magnetosphere being inside the stellar 
Alfven surface or at least in contact with it. If the planet is out- 
side the Alfven surface, the strength of the interaction is remark- 
ably reduced and the perturbations produced by the motion of 
the planetary magnetosphere cannot propagate back to the stel- 
lar surface because the stellar wind blows them away. Since the 
changes in the stellar field intensity and geometry along the ac- 
tivity cycle can modify the radius of the Alfven surface up to a 
factor of 1 .5 - 3, we expect a remarkable variation of the SPI sig- 
nal and even its disappearance during the time intervals in which 
the stellar Alfven surface retreats inside the planetary orbit. 
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